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Abstract 



We calculate the long distance contribution to D +, ° — > 7r + '°Z + Z~ decays by the use of 
a vector meson dominance model, in which the 0-meson plays the central role. The 
branching ratios obtained are 10 -6 and a few times 10 -7 for the resonance and non- 
resonance regions respectively. The analysis includes a calculation of D + — > n + (f), 
consistent with the experimental value. 



PACS number(s): 13.20.Fc 12.40.Vv 12.15.Lk 



1 INTRODUCTION AND THE MODEL 



New experimental limits for transitions involving change of flavour in charm meson 
decays were obtained recently at FERMILABJI], and at the Cornell Electron Storage 
Ring using the CLEO II detector 0. In particular, upper limits in the range 10~ 4 — 
were established for exclusive channels of the type D +, ° — > X + '°l + l~ , where X +, ° is a 
pseudoscalar or vector meson. 

The short distance process c — ► uj, driven by the magnetic penguin, is known to 
be of little significance for radiative decays of charm mesons as it leads to a branching 
ratio of 10~ 12 — 10~ n only, despite its enhancement by gluonic correct ions [Q. In decays 
to lepton pairs, the short distance c — > ul + l~ transition may be of more relevance, as 
it contains contributions from both form factors of the electromagnetic penguin, as 
well as contributions from the Z°-penguin and from the VF-box diagram. Indeed, 
the rate for the short distance c — > ul + l~ transition has been calculated []5| to be 
1.1 x 10~ 20 GeV, leading to a branching ratio of 1.8 x 10~ 8 for the inclusive process. 
Accordingly, exclusive process like D — > 7r/ + /~, expected to be approximately 10% of 
the inclusive rate, will have a short distance contribution to the branching ratio of the 
order of 10~ 9 . On this basis, it has been argued || that decays like D +, ° — > tt +,0 1 + 1~, 
— ► K + l + l~ constitute "... a large "discovery" window: seeing this decay occur at 
a branching ratio above ~ 10~ 7 would be strong evidence for new physics" . 

It is therefore of obvious interest to obtain reliable estimates for the long distance 
contributions to these modes. As one suspects that these contributions are the 
dominant ones in the Standard Model, one should have good control of their estimation, 
in order to perform a meaningful search for new physics in these decays. 

In the present paper, we consider the long distance contribution to the helicity- 
unsuppressed decays D — ► irl + l~, which supposedly are the better suited channel for 
checking the nature of the flavour changing neutral transitions. We note that presently 
the best experimental upper bounds for the branching ratios are as follows: 1.8 x 10~ 5 
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for £>+ -> 6.6xl(T 5 for J D + -> 7r + e + e-§, and 1.8xl(T 4 for D° -> 7rW~0, 

4.5 x 1(T 5 for D° -> 7r°e + e-[|. 

Our estimate for the long distance contribution is based on a vector meson domi- 
nance mechanism, similar to the approach widely used for calculating the long distance 
contributions in b — > sl + l~ decay Q. We present the details of our calculation for the 
D + — > 7r + Z + /~ channel, for which the experimental input required by our model is 
available. 

The basic assumption of our model is that the main long distance contribution to 
D + — > n + l + l~ is given by the transition .D + — > 7r + (V) — ► 7r + / + / _ , where V is a qq 
vector meson state. It is known that|| 

BR(D + ^n + (p) = (6.1 ±0.6) x 1(T 3 , 

BR(D+ -> 7T+p) < 1.4 x 10- 3 , (1) 
-> tt+cj) < 7 x lO- 3 . 
Moreover, the branching ratios of (J) for decays to lepton pairs are|| BR(<f) — > e + e _ ) = 
(3.00±0.06) x 10- 4 , 5^(0 -> //+//-) = (2.48±0.34) x 10~ 4 , while for p and the similar 
branching ratios are nearly one order of magnitude lower. Thus, we may restrict our 
considerations to the role of the 0-meson only, the contribution of p and u constituting 
a relatively small correction. 

Another kind of possible long distance contributions comes from the W- annihilation 



(or W-exchange) diagram, which is found for example to be large in B — ► P7|T0 
and needs therefore to be discussed further. In the present case of D — >• 7re + e~, the 
mechanism for the W-annihilation (or W-exchange) diagram is through D — * tc-/* and 
7* — > e + e~, where the virtual photon comes from one of the four quark lines. Because 
naive (constituent) quark model does not work for the light pseudo-Goldstone meson 
71, we need to use a hadronic model to estimate these contributions. At the hadronic 
level, the process D — ► 7T7* corresponds to the electromagnetic transition D — > -D7* 
followed by the weak transition D — ► tt, or D — ► 7r followed by 7T — ► 717*. These 
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electromagnetic transitions are described by the electromagnetic formfactors of the D 
and the 7r mesons, which can be taken in the VMD model as dominated by vector 
mesons as ip, p and u. The contribution from ip is smaller than those from the other 
two since this ip is highly off-shell. The remaining contributions from p and u in the 
W-annihilation (or W-exchange) mechanism correspond to the full processes D —>■ irp 
and D — > ttuj. These effects have been included in our VMD model and found to be 
small. 

In the second section we present our approach to the calculation of D + — > 7r + and 
in the third section we calculate its contribution to the D — > modes. 



2 THE D + 7T + d TRANSITION 



We begin with a treatment of the process D + — > 7r + 0, which is the main contribution 
to D + — > 7r + / + / _ in our approach. The effective Hamiltonian for this transition is 

Heff = ^V: s V cs a 2 sY(l - 75)^(1 - 7s)c, (2) 

where 02 is QCD-coefficient which is taken as \a 2 \ = 0.55 ±0.1 from the overall fit 
to nonleptonic D decays We shall assume the usefulness of the factorization 



hypothesis |pT|, thus 



< n + ( j ) \n eff \D + >= ^v: s V cs a 2 < 0|s 7 M (l - 7 5 )s|0 >< ir + \u^(l - l5 )c\D + >, (3) 
and we define by 



< 0|$7" S |o >= m^7- (4) 

will be determined from the observed lepton pair decay (f) — > e + e~. 
The hadronic matrix element is parameterized by two independent form factors f + 
and /_ as 



< 7T + |W7 M (1 - lb)c\D + >= U(pD+P7r)tM + f-(pD-p*)p- (5) 

The /_ form factor does not contribute in the decay to a lepton pair; / + is related by 
isospin symmetry to the form factor of the semileptonic decay D + — > ti 1 + j?i. Accord- 
ingly, we turn to the latter for learning the /+ form factor. 



Using the heavy meson chiral perturbation theory, Wise|]T2[ has calculated the /+ 
form factor near the zero recoil point q 2 m = {pD—Pn)m — { m D+ —m^o) 2 in D + — > 7i°l + i?i 
to be: 

/+«) = - ^7 (1-0- 



/d :i+^ d ~7), ( 6 ) 



2^2/^ m„- + A 

A = 777, £)• — ?7l£). 

The analysis of various existing experiments gives the following bounds on the param- 
eters in Eq. |13|1 : 



f D < 0.31GeV, 

(7) 

\g\ < 0.63. 



On the other hand, two measurement s [14] of the formfactor f+(q ) in the semileptonic 



2, , /2\ 1 -?1H 



decay have been analyzed under the assumption of a monopole behaviour for it, 

and the data may be summarized as[|] 

/ + (0) = (1.0 ± 0.3) x /f(0) = (1.0 ± 0.3) x (0.75 ± 0.03), (9) 

which shows the consistency of Eqs. (§D - ©, within the accuracy obtained so far. A 
new recent measurement by E687 gave |/+(0)//^(0)| = 1.03 ± 0.16 ± 0.02|15|], which 
is consistent with the value cited in Eq. (9). 
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Using now, for exemple, the reasonable set of values /o = 0.20GeV and g = 0.6, as 
well as Eq. (g), we find T(D + — > 7r+0) = (3.7±1.4) x 10 _15 GeV, which is in remarkable 
agreement with the experimental average|J ex PT(D + — > 7r + 0) = (3.8±0.4) x 10~ 15 GeV. 



This approach leads to /+(0) = 0.92 which agrees with the results of Refs. [[14], 15] . 

In the calculation of the leptonic decay D + — > 7r + / + /~ we shall normalize to the 
observed D + — > 7r + rate. Nevertheless, it is appropriate to emphasize at this point 
that the theoretical treatment described here accounts well for this mode; this provides 
the needed confidence in its use as the major contribution to the D — > 7r/ + /~ decays. 

3 LONG DISTANCE CONTRIBUTION TO D + 

TT+l+l- 

The effective coupling between </> and , via a photon propagator is 



90 = [e^- 7 > + ] (p+ ^_ )2 [^e^e 
= -=-ui-7 v i+7 , \2 e <p- 



(10) 



For the 0-meson decaying onshell, we replace l/(p+ + P-) 2 — * I/ 771 !- 

We shall assume that g^q 2 ) defined in Eq. (f|) does not vary appreciably with 



q 2 in the region of interest for our calculation ||16||, which is taken as m v < Jm\- < 
(rriD ~ m n), and we use for it the value determined from <fi — > e + e~ decay 

00 = (492MeV) 2 . (11) 
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The amplitude for the decay is then 

A{D+ -> ir+l+l-) 

--CL2V* S V CS < vr + |M7 M (l - 75 )c|D + > 



Gp 

9nv~ {p++P-)^P+ + P-)v/mle 2 _ „ (12) 



(p + +P-) 2 - m| + il^m^ 3 (p++p_) 2 



where we used for the propagator a Breit-Wigner form to account for the behaviour 
throughout the region of decay. 

Since at the 77-mass the sizable D + — > 7r + (77) — > ir + fi + fi~ opens [ BR(D + — ► tt + i]) = 
(7.5 ± 2.5) x 10~ 3 ], we shall impose a lower cut on the m + _ spectrum of the muon 
pairs above the r\ mass. We have also checked that changing this cut up to 700MeV 
does not affect practically our results. We shall make for convenience the same cut in 
the e + e~ channel. Then, the spectra for e + e~ and are essentially identical in the 

chosen region for m + _. By restricting our considerations to this region, we also avoid 
a possible ambiguity concerning th q 2 dependence in the region close to q 2 = 0. 

The contribution of <fi to the long distance D + — > 7r + Z + Z~ decay in the resonance 
region is given by 



dT(D + -> 7T+/+/-) 

r = / dm;_^ — 

'(m -ry2)2 am 
and we find it to contribute to the branching ratio 



— — 2 J -, (13) 



(Til) 



£#(£>+ -> 7T+/+/-) = 0.82 x 1(T 6 . (14) 



If we extend the limits in (13) to ± T^, this branching ratio becomes 1.22 x 10~ 6 . 

Now we turn to the region outside the resonance, which is the main object of the 
present investigation. We denote this contribution as T^r, and we calculate it from 
the amplitude of Eq. ([12] ) for the region m v < m H < (m^ — A) and (m^ + A) < 



m + _ < (mo — iTLn)- The branching ratios for the long distance contribution outside 
the resonance region, for several values of A, are given in Table 1. 



TABLE 1. The long distance contribution T^r to the D + — > 7r + fi + fi or D~ 
n + e + e~ decay rates, outside the ^-resonance (beyond ± A). 



A(in MeV) (in 10" 19 GeV) Branching ratio 

5 3.66 5.9 x KT 7 

10 1.92 3.1 x 10~ 7 

20 0.95 1.5 x 10~ 7 

40 0.45 0.73 x 10~ 7 



We remind the reader now that the short distance contribution || to these decays 
is about two orders of magnitude smaller than the values of Tables 1. Hence, the 
typical interference distribution in m + _ expected in b — > should not appear 

in D — > 7r/ + /~, the spectrum of the lepton pair in our case being given by the matrix 
element of Eq. ([12]) and phase space only. 

In the present paper, we concentrated on the D + — > 7r + / + /~ modes. For the parallel 
decays D° — > 7r°/ + Z~, we expect in our model a branching ratio smaller by approxi- 
mately 5, due to a factor of 2 from the AI = | weak D — > 7r transition and a factor of 
2.5 from the D + — D° lifetime difference. 

We remark that a previous long-distance calculation of D —>■ ttI + 1~ (first Ref. of 0) 
has considered it as evolving from the D — > D*ir process, with the virtual D* decaying 
to a lepton pair. A rate of 10~ 8 has been found for this contribution; this is consistent 
with our result as this specific diagram is obviously only part of the form factor we 
considered here for the transition (5). 



We conclude by stressing that our results, presented in Eq. (13) and in Table 
1, invalidate the earlier expectations || that branching ratios above ~ 10~ 7 would 
constitute a signal for physics beyond the Standard Model. 
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